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Abstract 

We shall proceed with the construction of normalizable Dirac wave packets for fermionic particles 
(neutrinos) with dynamics governed by a "modified" Dirac equation with a non-minimal coupling 
with an external magnetic field. We are not only interested on the analytic solutions of the 
"modified" Dirac wave equation but also on the construction of Dirac wave packets which can 
be used for describing the dynamics of some observable physical quantities which are relevant 
in the context of the quantum oscillation phenomena. To conclude, we discuss qualitatively the 
applicability of this formal construction in the treatment of chiral (and flavor) oscillations in the 
theoretical context of neutrino physics. 
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Since those early days when Dirac had derived the relativistic wave equation for a free 
propagating electron [1], several efforts have been produced in the literature to solve the 
Dirac equation with other analytical forms of interacting potentials, from central potential 
solutions [2, 3] to recent theoretical attempts to describe quark confinement [4, 5]. In fact, 
obtaining exact solutions of a generic class of Dirac wave equations becomes important since, 
for many times, the conceptual understanding of physics can only be brought about by such 
solutions. These solutions also correspond to valuable means for checking and improving 
models and numerical methods for solving complicated physical problems. In the context 
in which we intend to explore the Dirac formalism, we can quote the effect of neutrino spin 
flipping attributed to some dynamic external [6] interacting processes which come from the 
non-minimal coupling of a magnetic moment with an external electromagnetic field [7] and 
which was formerly supposed to be a relevant effect in the context of the solar-neutrino 
puzzle by suggesting an explanation for the LSND anomaly [8-10]. To be more specific, 
our aim in this manuscript is to try to accommodate the Dirac wave packet formalism 
[11, 12] and a further class of static characteristics of neutrinos, namely, the (electro)magnetic 
moment which appears in a Lagrangian with non-minimal coupling. By following this line of 
reasoning, we are not only interested in solving a "modified" Dirac wave equation but also 
in constructing Dirac wave packets which can be used for describing the dynamics of some 
observable physical quantities which are relevant in the context of the quantum oscillation 
phenomena. 

Despite their electric charge neutrality, neutrinos can interact with a photon through loop 
(radiative) diagrams [11, 13]. The Lagrangian for the interaction between a fermionic field 
and the electromagnetic field F^ u (x) is given by 

C = ^fi^(x)a fMU F^(x)^(x) + h.c., (1) 

with x = (t,x), a^y = | [7^, 7„] and F^ v (x) = d^A v {x) — d u A^{x) where we have not discrim- 
inated the flavor/mass mixing elements since we are initially interested just in solving the 
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equation of the motion described by 



i^(x) = {ol-p + [3 
at 



III ( ^//rr /(; , /••"",:,■, • ll.r. 



V(x) 

= {a- p + (3[m- Re{n)H ■ B{x) + Im{n)ot ■ E(x)]} V( x ), (2) 



where a = YH=i a ^ = 2~2l=i 7o7fc^> P = 7o> B(x) and -E(x) are respectively the magnetic 
and electric fields. The real (imaginary) part of /x represents the magnetic (electric) dipole 
moment of the mass eigenstate represented by ty(x). But it can be demonstrated [13] that, 
for Dirac neutrinos, the electric moment must vanish unless CP is violated, and, for Majo- 
rana neutrinos, the electric moment vanishes if CPT invariance is assumed. In the standard 
SU(2) L ® U(l)y electroweak theory [14], if a positive chirality neutrino is an SU(2) L sin- 
glet, the expression for fi can be found from Feynman diagrams for magnetic momentum 
corrections [13] and turns out to be proportional to the neutrino mass (matrix), 

3eG 3m e G O7v1fl _ l0 rn v 

a = — — — m = — = — u s m v = 2.7 x 10 a„ (3) 

8V2n 2 4V2tt 2 m N 

where G is the Fermi constant and uin is the nucleon mass. In principle, for m u w 1 eV, 
the magnetic moment introduced by this formula is exceedingly small to be detected or to 
affect astrophysical or physical processes. 

From the theoretical point of view, the physical implications of the non-minimal coupling 
with an external magnetic field can then be studied by means of the eigenvalue problem 
represented by the Hamiltonian equation 

H(p)V n (p) = {a-p + P[m-Re(n)?:-B]}V n (p), (4) 

which represents the time evolution operator of a spinor ty(x) for times subsequent to the 
creation (t — 0). The most general solution to the above eigenvalue problem is represented 
by 

E n (p) = ±E,(p) = ±\jrrf + p 2 + a 2 + (-l) s 2 ^m*a* + (p x af\ s = 1, 2, (5) 
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where a = Re{ji)B and the complete set of orthonormal eigenstates corresponds to 



where p a is the relativistic quadrimomentum, p s = (E s (p),p), N(p s ) is the normalization 
constant and 

Af = A 2 s (p) ± 2m\a\ - a 2 , af = 2E a (p)\a\ ± (A» + a 2 ) 

with 

A 2 s (p) = E 2 ( P )-(m 2 +p 2 ) + a 2 . (7) 

We can observe that the components of ^ ± {p 3 ) does not satisfy any auxiliary conditions, 
namely, at a given time t, they are independent functions of p and the eigenspinors does 
not represent neither chirality nor helicity eigenstates. Thus, if we seek the plane wave 
decomposition as 

^! + (x) = exp [—i(E s (p) t — p-x)] U(p a ), for positive frequencies and 

\P (x) = exp [ i(E s (p) t — p ■ x)] V(p s ), for negative frequencies (8) 

the time evolution of a plane wave packet solution ty(t, x) in this circumstances can be written 

as 

*(t,x) = f-^J2{b(p.)U(p.)exp[-iE.(p)t] 

+ d*(p s )V(p s ) exp [+i E s (p) t}} exp [ip-x], with p a = (E a , -p)(9) 

which, however, requires some extensive mathematical manipulations for explicitly con- 
structing the dynamics of an operator O as 



O(i) 



J M^OVM. (10) 



Meanwhile, if the quoted observables like the chirality 7 5 , the helicity S • p or even the 
spin projection onto the magnetic field X • B) commute with the Hamiltonian H, we could 



reconfigure the above solutions to simpler ones. To illustrate this point we shall limit our 
analysis to very restrictive spatial configurations of B so that, as a first attempt, we could 
calculate the observable expectation values which appear in Eq.(lO). Let us then assume that 
the magnetic field B is either orthogonal or parallel to the momentum p. For both of these 
cases the spinor eigenstates can now be decomposed into bi-spinors which are eigenstates of 
the spin projection operator cr • B as 



U( Ps ) = N + ( Ps ) 



(11) 



and 



V(p.) = N-(p s ) 



(12) 



<P (l>.) 
X'iPs) 

i. e. beside of being energy eigenstates, the general solutions U{p s ) and V{p s ) will become 
eigenstates of the operator £ • B and, equivalently, of £ • a. The Eq.(4) can thus be 
decomposed into a pair of coupled equations like 



(±E S - m + cr • a) <pf = ±<r-pxf, 
(±E s + m- a ■ a)xf = ±(r-p<pf, 



(13) 



where we have suppressed the p s dependence. The constraint imposed by the above equations 
leads to following representation for the energy eigenstates, 



U( Ps ) = N + ( Ps ) 



and 



V(p.) = N~(p.) 



(E s +m+(T-a)(T-p , 

\{E s +m)*+a 2 \ ^ s 



(E s +m+(T-a)(T-p _ 

l(E a +m)2+a 2 \ Xs 

x; 



(14) 



(15) 



Since </?+ 2 
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represent orthonormal eigenvectors of cr • a, we can 

J V 1 

immediately deduce the orthogonality relations 



uHps)U( Pr ) = vHps)V{ Pr ) = 5 sr , 

U\ Ps ) lo V( Pr ) = V\ Ps ) lo U{Pr) = 0. 



(16) 



The above results can be obtained by simply introducing the commuting relation [a ■ p, er ■ 
B] — which is derived when p x B = and the anti-commuting relation {a ■ p, cr ■ B} 
when p ■ B = 0. Now let us separately summarize other relevant information for both of 
these restrictive cases. If p ■ B = 0, the eigenspinor representation can be reduced to 



U(p s ) 



e + m 



2e 



and 



V(p.) 



m 



2e 



crp 4 

_ to+m *~P s . 



(TO 

— Y 

eo+m 



with e = a/p 2 + m 2 and the energy eigenvalues 



(17) 



(18) 



±E. = ±[e + (-iy\a\]. 



(19) 



The closure relations can be constructed in terms of 



- E V(p.) ® V^.JTo 



2e 
2e 



1 - (-l) s 7 oS • a 



E 



1 - (-l) s 7oS • a 



(20) 



where p = (e ,p). If p x B = 0, the eigenspinor representation can be reduced to 



U( Ps ) = 



E s + m s 



IE, 



and 



V(p.) 



with m s = m + (-l) s |a| and the energy eigenvalues 



(21) 



(22) 



±E S = ±v^T 



mi 



(23) 



In this case, the closure relations can be constructed in terms of 



U{ Ps )®U\ Ps ) lo = 



f 7^ + m s 



1,2 ^ 



1 - (-1) S E • a 



s = 1,2 



s = 1,2 



7 M p^ + m s 
2E« 



(24) 



Finally, the calculation of the expectation value of 0(t) is substantially simplified when 
we substitute the above closure relations into the wave packet expression of Eq. (9). To 
clarify this point, we suppose the initial condition over \I/(0, x) is the Fourier transform of 
the weight function 

if( P - Pi )w= i b (Ps)U(P*) + d*(ps)V(p s )} (25) 

s=l,2 

so that 

j2^<P(p-Pi)expbP-B\v> (26) 

where w is some fixed normalized spinor. By using the orthogonality properties of Eq. (16), 
we find [11] 

6(p.) = ip{p-p i )U ] { Ps )w, 
d*(p s ) = <p( P - Pi )VHp s )w. (27) 

These coefficients carry important physical information. For any initial state \I>(0, x) given 
by Eq. (26), the negative frequency solution coefficient d*(p 3 ) necessarily provides a non-null 
contribution to the time evolving wave packet. This obliges us to take the complete set of 
Dirac equation solutions to construct the wave packet. Only if we consider the initial spinor 
w being a positive energy (E s (p)) and momentum (p) eigenstate, the contribution due to 
the negative frequency solutions d*(p 3 ) will become null and we will have a simple expression 
for the time evolution of any physical observable. 

As we have already noticed, under the point of view of physical applicability, the above 
discussion can be useful in describing the dynamics of propagating neutrinos which non- 
minimally couple with an approximately constant external magnetic field B. In fact, the 
advent of the neutrino physics [16-18] has stirred up an increasing number of theoretical 
works where the quantum oscillation mechanisms [15, 19-22] have been deeply analyzed. 
In particular, the correspondence between helicity and chirality has frequently caused some 
confusion in the literature where the concept of helicity has been erroneously used in the 
place of chirality for particles with m = 0. In fact, some conceptual differences between the 
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chiral and helicity operator has been analyzed [23, 24] in the context of the dual-helicity 
eigenspinors of the charge conjugation operator (Elko) [26]. 

In the standard model of flavor-changing interactions, neutrinos with positive chirality are 
decoupled from neutrino absorbing charged weak currents [25]. Consequently, such positive 
chirality neutrinos become sterile with respect to weak interactions. By reporting to the 
formalism with Dirac wave packets [11, 12] which leads to the study of chiral oscillations 
[25] (in vacuum), with the formal procedures presented in this manuscript, one could obtain 
the immediate description of chiral oscillations and spin flipping in terms of the Hamiltonian 
(2) dynamics by recurring to the following dynamic equation 

|( 7 5 ) = 2i( 7o75 [m-iHu)£.B]> (28) 

and 

j f (h) = iiteM< 7 o(E x B)) ■ p, (29) 

where the particle helicity is defined as the projection of the spin angular momentum onto 
the vector momentum direction h = |S p. These results are obtained by means of a couple 
of very simple calculations, nevertheless, we believe that, in spite of its simplicity, it is im- 
portant in the large context of quantum oscillation phenomena. We still remark that, in the 
standard treatment of vacuum neutrino oscillations, the use of scalar mass-eigenstate wave 
packets made up exclusively of positive frequency plane-wave solutions is usually (implic- 
itly) assumed. Although the standard oscillation formula could give the correct result when 
properly interpreted, a satisfactory description of fermionic (spin one-half) particles requires 
the use of the Dirac equation as evolution equation for the mass-eigenstates. Consequently, 
the spinorial form and the interference between positive and negative frequency components 
of the mass-eigenstate wave packets lead to the possibility of chiral coupled with flavor os- 
cillations [15] which, eventually, deserves a very careful investigation for cases of neutrinos 
which propagate through external interacting (magnetic) fields. 

To summarize, we would have not been honest if, by considering an arbitrary spatial 
configuration for the magnetic field, we had ignored the complete analysis of the general 
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case comprised by Eqs. (4-6). Meanwhile, such a peculiar observation leads to the formal 
connection between quantum oscillation phenomena and a very different field. It concerns 
with the curious fact that the above complete (general) expressions for propagating wave 
packets do not satisfy the standard dispersion relations like E 2 = m 2 + p 2 excepting by the 
two particular cases where E s {p) 2 = m 2 s +p 2 for p x B = or e 2 = m 2 +p 2 for p ■ B = 0. By 
principle, it could represent an inconvenient obstacle which forbids the extension of these 
restrictive cases to a general one. However, we believe that it can also represent a starting 
point in discussing dispersion relations which may be incorporated into frameworks encoding 
the breakdown (or the validity) of Lorentz invariance. 
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